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The Reynolds-stress and triple product Lag models were created with a normal stress 
distribution which was defined by the accepted 4:3:2 distribution of streamwise, spanwise 
and wall normal stresses, and a ratio of 7, = 0.3k in the log layer region of high Reynolds 
number flat plate flow, which implies R}, = wars x 2.96. More recent measurements show 
a more complex picture of the log layer region at high Reynolds numbers. The first cut 
at improving these models along with the direction for future refinements is described. 
Comparison with recent high Reynolds number data shows areas where further work is 
needed, but also shows inclusion of the modeled turbulent transport terms improve the 
prediction where they influence the solution. Additional work is needed to develop a 
model that better matches experiments, but there is significant improvement in many of 
the details of the log layer behavior. 


Nomenclature 


free stream Mach number 
Reynolds-stress tensor [= Uu;U 


- - _. al Ou; Ou; 
strain rate tensor [= 5 ( Dee, + Bx, )| 


Turbulent transport tensor [= U;U Up 


rotation rate tensor [= 4 (3 _ a )| 
a 

molecular viscosity 

eddy-viscosity 


kinematic viscosity [= u| 


specific dissipation rate [= a 


Favre average Cartesian velocity fluctuation components 
mass density 
wall shear 


. . . _ ‘i , - 
turbulent kinetic energy per unit mass [= sujet] 


homogeneous turbulence dissipation per unit mass [= 8*kw] 
static pressure 

free stream velocity 

instantaneous Cartesian velocity components 

Cartesian position coordinates 

Favre average Cartesian mean velocity components 
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I. Introduction 


Accurate computational flowfield predictions are essential for both design and operation of aerospace 
vehicles. As computer speeds and memory size continue to increase, Computational Fluid Dynamics (CFD) 
can be used to predict the flowfield around not only simple shapes but also complete vehicle configurations. 
The advances in computer clock speed and memory capacity have allowed the modeling of turbulent flow, at 
lower Reynolds numbers, using Direct Numerical Simulations (DNS). Large Eddy Simulation (LES) continues 
to be developed for application at higher Reynolds numbers, but for complex configurations, DNS or even 
LES are still impractical because the grid resolution required (in both time and space) is well beyond current 
computational capabilities. 

Reynolds-stress turbulence models were envisioned to overcome a number of shortcomings evident in 
simple Boussinesq eddy-viscosity models. Although Reynolds-stress models have had a long history of devel- 
opment,!’? they have had, until recently, limited success in actually overcoming these limitations in practice. 
Reynolds-stress models have enjoyed a resurgence in the past few years,” ° with one new methodology incor- 
porating the desired flow history effects on the Reynolds-stress tensor in a formulation that is numerically 
robust.°'’ This Lag methodology®*!? allowed a further expansion of the flow history to include triple ve- 
locity products in a bid to obtain more accurate and complete turbulent transport predictions. This new 
model,!° denoted “TTR” for Turbulent TRansport, augments the second-moment predictions of the Lag 
Reynolds-stress models, adding field equations for the third-order-moments. These are an attempt to fulfill 
the need for turbulent transport predictions in regions of separation, where their relative importance is larger 
than it is for attached flows. 

Previous work on these models was focused on fleshing out the mathematical system that would 
allow reliable, robust computation of complex flowfields with the full benefit of Reynolds-stress and higher 
moment models. The stress-strain relationship utilized in that effort was based upon a classic, but inaccurate, 
outdated set of relations: 


6, 10-12 


Ty = 0.3k 
(Ri: R22: R33) = (4:3:2) 


There has been much work on more complex relations in the past few decades,!*!® and a final consensus is 


yet to crystallize, but there is little doubt that, for instance, the streamwise Reynolds stresses (Rip were 
at least double what this relationship would predict. To be clear, the Lag models are not alone in retaining 
this out of date turbulent paradigm,” but rather appear to be the first to attempt to embody the newer 
experimental knowledge. 

For this paper, the main changes have to do with flat plate skin friction predictions (which are now 
believed to be 6% too high!” !°) and law of the wall constants (from « = 0.41, A = 5 to Kk = 0.385, A = 4.2) 
but more importantly the rea for a flat plate boundary layer. The earlier dogma yields an eel = 2.96, 
where the more modern data!+?°?? show Rf, > 5.5 . Initial investigation of the effects of raising this 
ratio showed that it did affect separation predictions slightly, but had a substantial effect on reattachment 
predictions for the Bachalo-Johnson bump, causing an earlier reattachment more in line with experimental 
data. Comparison with relatively recent experimental data,?? which was kindly provided by the first author 
of that report, suggests that more work is going to be required to obtain a Reynolds-stress model that will 
predict the actual log law behavior of the Reynolds stresses, so separation prediction comparisons with this 
interim version would be premature, even though they are promising. 

As the triple products depend directly on the Reynolds-stress field, increasing R,; levels will generally 
increase T;;; levels, and as the triple products appear to give a more accurate prediction of the turbulent 
state in regions with low production such as separated zones, the triple product models are a central focus 
of this paper. The low Reynolds number near-wall, region, yt < 100, is not a concern of this effort. The 
Wilcox-1988 k — w model has been tuned to give good predictions of attached wall shear layers for both 
favorable and adverse pressure gradients, and this is retained as the “near-wall” model. More complex 
models which would handle the physics extant in this region are envisioned and under development, but are 
not part of this paper. 

The Reynolds-stress and triple product models described in this paper are not the final product, and miss 
a number of key features which were revealed with the more recent high Reynolds number data. However, 
they are closer to that data than the previous versions. This paper describes the predictions of the TTR 
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model in this chrysalis version. In line with the improved predictions of the TTR model on the rotating 
pipe flow,!? the triple product terms do seem to improve the model predictions in the regions where they 
are active, in this flowfield at the laminar/turbulent interface of the flat plate flowfield. 


II. Experiment Description/Computational Methodology 


II.A. Experiments 


Detailed experimental results of the Reynolds stress R}, are taken from two sources. The first source 
is the Driver CSO and BSO flowfields,?* where the inflow boundary layers 0.3m upstream of the spinning 
cylinder and adverse pressure gradient have been used. This flowfield is well documented, with skin friction 
determined by oil flow interferometry and three component LDV velocity measurements. The CSO flowfield 
is an axisymmetric flowfield, but the boundary layer height to cylinder radius is small, and it is treated here 
as simply a source of high quality flat plate data. 

The second source is from Marusic,?” and is the tripped boundary layer data (SP40) discussed in that 
paper. This is another well documented flow field, in this case the floor of a wind tunnel. The skin friction 
(or more correctly u,) was determined by a fit to a composite velocity profile(one created to match the 
more recent and extensive flat plate data), and streamwise velocities were measured with hot wires. The 
Reg range of this dataset is extremely high compared to most experiments of the previous century, With 
the exception of a few experiments!* !”!*:?4 which were completed at its conclusion. As high a Reynolds 
number as these experiments reach, it is noted that it corresponds to the Reynolds number at the first 
passenger window on a business jet. Measurements extend deep into the boundary layer, into the sublayer 
and production peak region. For TTR model development, the details for yt > 100 are the object of 
interest. The freestream turbulence level reported for this tunnel will be used with the “Wind Tunnel” 
boundary conditions described below. The turbulence intensity for this tunnel is reported as wu’ /Uoo < 0.002, 
and the measurements of u//u, at the boundary layer edge are consistent with this. A selection of Ry, 
measurements from other experiments!*!%!* are also compared to the model predictions. 


II.B. Computation 
IT.B.1. Grid 


The flat plate grid used 
(Fig. 1) is 513 x 513, with 
initial wall normal spacing 
of 10-82 and an initial ax- 
ial spacing of 10-4Z. The 
wall normal stretching ratio oe 
is less that 1.03, and the ax- eee 
ial stretching ratio is less that Weeeeie tees eeceeee rece 

1.02. The grid vertical extent ae see cee oeee: 
starts at the leading edge at em oS goof geese see 

0.01L, and linearly grows to 
0.3L at the trailing edge, very 

similar to grids from.'° The 

grid is sufficiently fine to allow Figure 1: Flat plate grid; blue lines are every 4th grid point. 
simulation of this flowfield at 

axial locations that are less than 0.005L, which corresponds to a Reynolds number based on run distance 
and freestream velocity of 500 x 10°, nominally the lowest Reynolds number at which it would be possible 
to have turbulent flow at normal conditions. 


IT.B.2. Boundary Value Problem Definition 


The flowfield to be studied is the canonical low Mach number turbulent flat plate. This simulation is 
accomplished with a fully turbulent plate with a Reynolds number based on plate length of 100 x 10°, which 
gives a long region of fully developed turbulent flow—certainly 90% of the plate length. The simulation 
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Mach number was set to 0.2, yielding an essentially incompressible flow field without requiring low Mach 
preconditioning. 

The boundary conditions on the plate are viscous, adiabatic wall along z = 0, characteristic boundary 
conditions along the inflow plane and upper edge, and simple extrapolation along the exit plane. A detail 
that is generally not discussed at any length is the freestream conditions for k and Rr = rae The standard 
conditions used are to set kj. = 1 x 107°, and Rer = 0.1, which is essentially laminar freestream. The 
turbulent kinetic energy continues to decay from the inflow edge, and k and Rev are actually lower for any 
part of the solution domain without significant shear strain. This can be thought of as “Flight” conditions, 
where the atmospheric turbulence is vanishingly small, and this is the usual boundary condition imposed. 

An alternative is to attempt to match the turbulence state existing in a wind tunnel test section, where 
the grid turbulence of the last screens in the settling chamber is accelerated through the contraction section, 
then traverses down the tunnel at a nearly constant freestream turbulence intensity. This is done by setting 
k to the desired level, and Rer > 1000. If a measurement of the decay of k down the test section is available, 
this can be used to fix Rey, but in this paper it is set to 1000 for the cases where non-zero freestream 
turbulence is being simulated. This value gives a nearly constant freestream k over the length of the plate. 
Choosing higher levels of Rep give results that are essentially similar. 

For most of the paper, the standard “Flight” turbulence level boundary values are used, but some results 
look at the effects of freestream turbulence on the model, and for those the “Wind Tunnel” method is used, 


with “|. = ,/2#|.. noted. 


II.B.8. Numerical Method 


The code used in this study was a modified version of OVERFLOW 2.2k.?°:?° The modifications included 
the addition of Lag, Lag-R;;, and TTR models along with the high speed modifications.?’ Matrix dissipation 
was used with smoothing parameters as recommended by earlier studies of high-speed flows with this code?’ 
with one critical change. The critical difference in the matrix dissipation smoothing parameters used is that 
the eigenvalue limiters are set to zero. Matrix dissipation?® is appropriate for this flowfield. 

Results from the HLLC scheme, as coded in OVERFLOW,?””:*° were compared with the matrix dissipation 
scheme results, and the results agree with the modified matrix dissipation scheme. The Pulliam-Chaussee 
diagonal scheme,*! with variable time stepping or a constant Courant number (CFL) and multigrid was 
used as the relaxation method. Grid sequencing (called full multigrid in OVERFLOW) was utilized, and 
allowed a check on grid convergence as well as drastically reducing the CPU time required to fully converge 
the results 

In general, spatial convective terms and diffusion terms were all second-order accurate. For the modeling 
of the convection terms of the turbulence models, second-order upwind was used on all the Reynolds Averaged 
Navier Stokes(RANS) models. The Lag methodology does require second-order accuracy (or better) since the 
field equations defining the lagged turbulent variables are a balance of convection and source with no diffusion 
terms by design—purely hyperbolic equations to accurately mimic the history effects so clearly evident in 
turbulent flow. In the turbulent transport level equations, not all the equations are purely hyperbolic, but 
they are all more driven by convection terms than standard one or two equation models. 


ILB.4. Turbulence Models 


For this paper, two Lag-T;,;, Triple product/Reynolds stress models are investigated. The baseline model is 
the one described in earlier papers, and includes the modifications found necessary to simulate the rotating 
pipe case.!? On the flowfields described in this paper, the results obtained are essentially the same as those of 
the corresponding Reynolds stress model (“926”).” The second model described here is based on that model, 
but with adjusted constants to obtain substantially higher A values in the flat plate log layer to better 
match the experimental results!* 2° ??:°?:°3 which show the much more complex behavior of the Reynolds 
stresses in flat plate boundary layers. Not all of this complexity is captured in this model. Regardless, 
the predictions are closer to experiment than the un-modified (baseline) model, and these revisions are a 
stepping stone to a model that captures the more complex behavior more completely. 

In tuning these constants, several constraints were maintained on the turbulence model. One of the more 
critical, and often overlooked, is that the model give the correct decay rate for isotropic turbulence.” This 
particular constraint is important for any model that is attempting to calculate separated flowfields. In the 
separation region, the turbulent production becomes small, and having the decay rate correct ensures that 
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the dissipation, which is no longer matched by production, has the correct physical behavior. This constraint 
was maintained in the models described here. 

The TTR model is built on top of the preceding Lag methodology models, and the equilibrium two 
equation turbulence model that is at their heart, the Wilcox 1988 k-w model. The third-order-moment 
model comes from the exact Reynolds-stress and turbulent transport equations. 


Or (Rij) + On (ue Riz) = — RjrOnUi — RixOgUj — OnTijn + VOKOn Ri; (1) 
+13 — 270; (uj) Ox (u5) 

On (Tijn) + Or (mijn) = — Titi — TinUs — Tein; (2) 
+ RAR +R RARGT + RAR; 
+ VOAT jk 


+Mligk — O1(Qijat) — €ijx 
where the neglected red terms are 


1 
Thay = [wi + i101) (3) 
1 
Thigh 5 [un + wie OW) + RH O;7)) (4) 
Qijet =UZU; Uj, UY (5) 
Eijk =2v (Wa Au) + ui O1(uy )O(ujy) + uh Oi(w, )Ou(ul)) (6) 


The TTR turbulence model including turbulent transport (Zj;,) terms is: 


O; (pk) + O; (ourk) =p [Riz Siz — B* kw] +O; ((u + oer) Ak) — Agd; (pTiit) (7) 
A: (pw) + (puw) =apS* — Bow? + A ((u + oop) Ow) (8) 
Ot (pRiz) + O1 (ow Rij) =Aopw (a = Ris) (9) 
Ot: (pTign) + O; (puTijn) =Aopw (ae = Tijn) (10) 
where 
e A — = _ 
aa = aa [TitiTn + TAU; + Trin; — Riz ARw — RjeARU — ReiARyi] 
1 
+ Anite [Or («+ Oren) OTijx)| (11) 
Mty = pk/w V 2 
S = 4/2 (SigSijy — Six /3) 
P = RiSis 
7 : 1 Ou; Ou; 
i as Sty = 5 Ox; = Ox; 


Most of the parameters for this model are set by the requirement to retain the equilibrium predictions of 
the underlying k-w model. The equilibrium Reynolds-stress relation is one of the three described in earlier 
Reynolds-stress model work,’ denoted as the “926(Redistribution)” Equilibrium Reynolds-stress relation. 
This constitutive relation is most directly related to the explicit algebraic Reynolds-stress models.°** °° The 
terminology is borrowed from the paper introducing this relation,’ and contains production terms which are 
not the full Reynolds-stress production terms, but do yield log layer anisotropies consistent with the classic 
4:3:2°" relation. These “production” terms (which are actually only redistribution terms, and none of the 
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work terms of actual production) are 


Pu = 2(k941/Ai + (Raz 231 — Rig 212)) (12) 
Po. = 2(kS22/A1+ (Rig Q12 — Reg Qe3)) (13) 
P33 = 2 (k933/A1 + (Rez Qe3 — Rei 231)) (14) 
with corresponding off-diagonal terms 
Pig = 2kS12/A1 + (Reg 31 — Rei 2o3 + (Riz — Rag) 212) (15) 
Po3 = 2kS93/A1 + (Rei Qig — R1g231 + (Rog — R33) 29) (16) 
P31 = 2kS31/A1 + (Rig Qe3 — Rog Q12 + (R33 — R11) 231) (17) 


then the equilibrium Reynolds stress is given by: 


a A Ls 
Raj) = Zhbig — — (Pay — 3 PSs) 
1 
of Aas (Or ((u + opr) ORi;) — As% (pTi;1)| (18) 
0 pw 


The model parameters which are constant for all TTR models in this paper are: 


Ag = 1.0 A, = 1.0 A3 = 1.0 A, = 1.0 
On = 0.3 or = 1.0 o, = 1.0 


For the model denoted TTR(original), the remaining parameters used are: 
Aj == a = 0.35 g= 04 oS U12 63 = 0.3 


while the model denoted TTR(revised) uses the parameter set 


A, = fe a = 0.008 B = 0.00375 B* = 0.0045 Ow = 0.3166322 
Note that even though 6 and §* have been altered from the 1988 k-w choices, their ratio is retained providing 
the same isotropic turbulence decay rate in all models. 

The parameter Aj essentially chooses R7, in the log layer. Increasing A, increases the R}, log layer value. 
The remaining parameters that are adjusted are chosen to maintain the isotropic decay rate discussed above, 
to match the log law velocity profile for a flat plate (there is a relationship between these parameters which 
fixes «, much like in the analysis for the underlying k -w(88) model’), and to match skin friction behavior 
expected for flat plate flow fields. There remains one free parameter which can be adjusted to match the 
flat plate skin friction relation. 

The constant set (in red, which matches the color used in Fig 2) with A; = 3 is labelled as “T'TR” in the 
figures, and is essentially the original TTR model, with some minor constant adjustments that were required 
to work in the rotating pipe flowfield.‘° So far, these slight changes have not introduced any significant 
differences in other flowfields, including the one reported in this paper. 

The constant set (in green, again matching the color in Fig 2) with A, = 4 is labelled as “TTR(Revised)” 
in this paper, and while this is a stepping stone on the way, this is not expected to be the final version of 
the complete model. Significant improvement in details of the turbulent flowfields in wall bounded flows has 
been made, but there is still work yet to be done to obtain a model which matches the current experimental 
knowledge, which has expanded and improved greatly in the past few decades. 

A full exploration focusing on separated flowfields is planned, and that paper will include a full comparison 
on a wide selection of flowfields, with an emphasis on separation and reattachment predictions. That work 
will be the model which will be used to meet the goals of the 40% improvement milestone of the Revolutionary 
Computational Aerosciences program.°® 
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III. Results 

The flat plate is the primary canonical wall- 
bounded turbulent flow field. Any prediction 0.005 [S717] TTT a eS 
method that deals with turbulent flow over walls — TTR 
should be able to give reasonable predictions of this . ———- TT R(Revised) 
flow field to have any actual utility. For eddy viscos- ja i) 46% — Lag-up 
ity models, this comes down to providing both skin a —— Karman-Schoener 
friction predictions as well as matching the axial ve- ne Nagib 
locity profile, the law of the wall. . © Exp - Marusic(2015) 


Recent experimental work has slightly altered 0.004 
both the skin friction and the « value at the heart 
of the law of the wall. In previous work,°:>!? the 
Karman-Schoenherr correlation of skin friction to 9 9935 
Reg was used. It was understood that there was 
an uncertainty in the experimental data of at least 
2%,°° but more recent work!’ has was able to ob- 
tain a much better agreement when accurate, inde- 
pendent measurement of the skin friction was ob- 
tained. Figure 2 shows both the original Karman- 
Schoenherr correlation, and Nagib’s modified ver- 0.0025 
sion of Coles-Fernholz(designated “Nagib” in that 
figure): cy(Reg) = 2 (In(Ree)/K + C)~”, with « = 
0.384,C = 4.127. The revised model can be ad- 
justed to give skin friction reasonably close to the = 
new correlation, at least up the highest flat plate Littl poiiitl poitit 
Reynolds number data available in a ground based 103 104 105 
facility, shown by the rightmost symbol on this plot. 
The revised model’s cy vs. Reg predictions would 
get further and further away from the new correla- 
tion, assuming that it describes the data at higher and higher Reynolds number. For Reynolds numbers that 
would be involved in aircraft wing design, it is probably acceptable, but not for simulations of the entire 
aircraft, especially those including the fuselage. Another issue is the behavior at low Reg, where the model 
will not sustain turbulence for Reg < 1500, a point we will return to later. For now, recall that results shown 
in this figure were obtained with “Flight” (vanishingly small freestream turbulence) conditions . 

The symbols plotted in Fig. 2 are derived from data obtained from,?? and the mean velocity and Bay 
profiles associated with them are shown in Fig. 3. Also included in this plot are two inflow condition profiles 
from the Driver CS0/BSO0 flowfields. The lower portion of this plot is the velocity profile in inner coordinates 
(yt, ut). In all previous work, the Lag models have been tuned to fit the Coles log-law constants, « = 0.41, 
C =5. Newer data from two separate facilities with both oil flow interferometry and MEMS instruments, !” !® 
along with velocity measurements give convincing evidence*”’ that « = 0.385 + .005,C = 4.17. The revised 
model fits the new log law constants. The revised model does do a better job of fitting the wake region, for 
both the lower and higher Reynolds number data. 

The bigger difference between the baseline and revised model is shown in the top plots of Fig. 3, the 
normal axial Reynolds stress in inner variables, R{, vs yt. The baseline TTR model (and the Reynolds 
stress models with which it is associated) have log layer Reynolds stresses that are less than half of those of 
the revised version. Even with this dramatic increase of Rj, in the log layer, it is still lower than that seen 
experimentally in the lowest Reynolds number data that has a wide enough log layer to produce a “shelf” 
in R],-red o—data that corresponds to Reg = 17.2K. This is the region at the start of the log layer, where 
0; Ri, 0. For the highest Reynolds number experimental data, Reg = 35 x 10°, this is 70 < yt < 500. 

The details for y* < 100 are not expected to be matched, it is the log law region yt > 100 that we are 
attempting to model. All the details below yt < 100 are subsumed by the underlying k — w(88) model, 
which gives good skin friction predictions for favorable and adverse pressure gradients in attached boundary 
layers. This is a low Reynolds number region, the domain of sweep/ejection events, apparently large coherent 
structures, and other details that will be dealt with at later. Figure 4a focuses on the region of interest, 


0.003 


0.002 


Figure 2: Flat Plate Skin Friction 
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Figure 3: Velocity(lower plot) and Axial Reynolds-stress(upper plot) in inner variables — 7,, match. 


and shows some of the experimental data available from the more recent experiments. This does support 
the existence of a shelf, but the shelf level R}, as a function of Reg is not obvious. Note that a 6% error in 
estimating 7,, will lead to a 6% change in R},, and the recent changes in c estimates” are of this magnitude. 
Care was taken to use data that was believed to be scale resolved in this figure - indeed the data from Hites!” 
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Figure 4: Axial Reynolds-stress in inner variables(focus above sublayer) — Reg matching. Experimental data 
from M — Marusic,?? H — Hites,!” O — Osterlund,!® D —DeGraff'* 


is the data taken with two different hot wire sensors which are certainly self consistent. 

The other striking miss in prediction for Rj, is that the width of the log law region, and the boundary 
layer itself, are too large in this revised model. Surprisingly, the revision actually is an improvement over the 
baseline model, as can be seen in the improved predictions at the boundary layer edge. The triple product 
terms in these models can seem, in some cases, to have an anti-diffusive character, and as the magnitude 
of these terms is proportional to the magnitude of the Reynolds-stresses, it is consistent that the higher 
Reynolds-stress levels of the revised model could be providing more non-diffusive assistance. One detail 
that bears mentioning is that it does not matter whether 7,, or Reg is matched between experiment and 
computation. Figure 3 matches 7,,, where Fig. 4a attempts a match with Reg. As might be expected since 
the model has been adjusted to match the variation of cf with Reg, the conclusions about the model’s Ry, 
shelf extent being wider than experiment is evident in both figures. 

A detail that the model is getting right is the decay from the log region plateau to the freestream floor 
in ae. This can be seen in the upper right plot in Fig. 3. The plateaus are too wide, granted, but the 
slope of the right edge of the model prediction for cf = 0.00232 red curve lines up beautifully with the right 
edge of the experimental data for cr = 0.00218, and the experimental slopes (0,+ Rj) at the three different 
axial locations are very similar as the freestream is approached . Note further that this is not the case 
for the baseline models in the upper left plot of Fig. 3, in this plot, the slope of the curves do not match 
as freestream is approached. This slope is directly affected by the triple product terms of the model. In 
the flat plate flowfield, the regions where the triple product terms have any effect are at the edges of the 
turbulent region, in the sublayer, and in this region at the laminar-turbulent interface at the outer edge of 
the boundary layer. 

One of the more intriguing developments occurred as we attempted to tune the model to even higher 
Ry},. The model was able to reproduce the log law velocity profile and follow the cy vs. Reg curve, but 
the turbulence was not able to sustain itself until Re, > 2 x 10°, a result that would be impressive in a 
wind tunnel used for transition research. However, up to this point, all of these simulations had been run 
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with free stream turbulence levels that are more like flight than wind tunnels, that is they are extremely 
low turbulence intensity. As a final check, this version was run with a freestream turbulence level that was 
consistent with a high quality (though not transition research grade) wind tunnel - a turbulence intensity 
level of about 0.2%. Earlier modeling work suggested that this would not significantly affect the point at 
which turbulence flow could be sustained by the model, but that is exactly what occurred. 

With this surprising result, simulations of 
the TTR(Revised) model were conducted with 
freestream turbulence conditions that were more con- 
sistent with wind tunnel test sections, the alternative oN, 
boundary condition values described above. Figure 5 ce . / 0.005 
shows skin friction vs. run length Reynolds number, 
and behavior for plausible freestream turbulence in- 
tensities is consistent with what one would expect 
with wind tunnels of varying freestream turbulence : 
intensities. In this plot, the green curve (“Flight”) is 0.003 |- * + 
exactly the same simulation as is plotted in Fig. 2, E \ 
except that the skin friction is plotted as a function \ 
of Re, rather than Reg. Increasing the freestream \ 
turbulence level significantly lowered the run length \ 
Reynolds number Re, at which turbulent flow could : 
be maintained. These results were insensitive to the 9-9025 A 7 
value of Rer for Rer > 1000, as these all provided 
a roughly constant level of & over the length of the S 
plate. The “transition” Reynolds number thus ob- N 
tained is consistent with what might be expected with N 
the turbulence intensity changes in a wind tunnel. 
The TTR(Revised) model formulation is not unique —°°°°? 
in this property,*! and in hindsight, it might have 
been anticipated, though these two modeling formu- 
lations are certainly distinct. 

Having the model depend at low Reynolds num- Figure 5: Effect of Freestream Turbulence level on 
ber on the freestream conditions is a mixed blessing. TTR(Revised) model skin friction predictions 
It is not inconsistent with the model starting to pick 
up more and more details of the flowfield, but it adds a complexity to model development, another free 
parameter that needs to be assessed and understood. It may be necessary as more and more of the physical 
processes of the turbulent boundary layer are modelled to include this detail. 


~ 0.002 


IV. Conclusions - Further Work 


The modifications to the Lag-T;;, (TTR) model are continuing. While the model does a better job of 
matching recent high Reynolds number data, there is work yet to be done in getting some of the details that 
are emerging. The existence of a plateau in By (and presumably the other R;; in the log layer) is a welcome 
development. If the Reynolds-stress distribution in the log layer had a more complex behavior, matching 
that behavior would be much more difficult. 

The current model versions miss the data in two major ways. First, the model in its current form predicts 
a single Rf, at sufficiently high Reynolds number, and experimental data suggests that the plateau levels 
increase with increasing Reynolds number. Second, the plateau extent of the current formulation is overly 
broad, extending too high in yt. The plateau size predicted by the model does increase with increasing 
Reynolds number, and with the triple product terms, it does have the correct decay in ce as the freestream 
is approached, so these details will need to be retained in the further development. A final point to be 
stressed is that the experimental picture (e.g. rea plateau levels as a function of Reg) is still not completely 
clear, and the model development is completely dependent on accurate experimental data. 

The transition behavior is intriguing, but adds additional details to check, as might be expected as 
more physical processes are added to the turbulence model. In a flat plate flowfield, the sublayer and 
laminar/turbulent interface regions are the ones most affected by the addition of the triple product terms. 
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The TTR model revisions improve in prediction of the laminar/turbulent interface region. There are two 
areas in which the current model need additional improvement. The width of the plateau region of Ry, 
needs to be reduced, and the R7}, plateau level variation needs to match that of experiment. Regardless, the 
general improvement in prediction of the Reynolds-stress terms does seem to translate into an improvement 
in T;;~ term behaviour, based on the laminar/turbulent interface region. 

A relatively famous statement was put forward by Joe Marvin*? thirty five years ago: “No single tur- 
bulence model emerges today that applies generally to the variety of flows encountered in Computational 
Aerodynamics: thus turbulence modeling remains an important pacing item”. Progress has been made, 
but this is still an accurate assessment today. The most important pacing item for turbulence modeling is 
detailed and complete experimental data from flowfields that can be accurately and completely modeled in 
CFD. The recent measurements at extremely high Reynolds number give a more detailed of the log layer, 
and one that suggests that there may be a field equation RANS model that can be devised to match that 
behavior. 
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